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I. INTRODUCTION 



The discussion of the electrodynamics of moving media has a long history. At present, the 
general structure of classical electrodynamics appears to be well established. Inparticular, in 



rol, the electric 



the generally covariant pre-metric approach to electrodynamics 
charge and the magnetic flux conservation laws manifest themselves in the Maxwell equations 
for the excitation H = (P, 7i) and the field strength F = {E, B), namely dH = J, dF = 0. 
These equations should be supplemented by a constitutive law H = H{F). The latter rela- 
tion contains the crucial information about the underlying physical continuum (i.e., space- 
time and/or material medium), in particular, about the spacetime metric. Mathematically, 
this constitutive law arises either from a suitable phenomenological theory of a medium or 
from the electromagnetic field Lagrangian. It can be a nonlinear or even nonlocal relation 
between the electromagnetic excitation and the field strength. The constitutive law is called 
a spacetime relation if it applies to spacetime ( "the vacuum" ) itself. 

Among many physical applications of classical electrodynamics, the problem of the inter- 
action of the electromagnetic field with matter occupies a central position. The fundamental 
question, which arises in this context, is about the definition of the energy and momentum 
in the possibly moving medium. The discussion of the energy-momentum tensor in macro- 
scopic electrodynamics is quite old. The beginning of this dispute dates back to Minkowski 
[3], Einstein and Laub [sj], and Abraham 9]. Nevertheless, up to now the question was 
not settled and there is an on-going exchange of confiicting opinions concerning the validity 
of the Minkowski versus the Abraham energy-momentum tensor, see, e.g., the review 10|. 
Even experiments were not quite able to make a definite and decisive choice of electromag- 
netic energy and momentum in material media. A consistent solution of this problem has 
been proposed in 4, lUl (cf. also the earlier work 12|) in the context of a new axiomatic 
approach to electrodynamics. 

Recently Feigel ISj has studied, theoretically, the dynamics of a dielectric magnetoelectric 
medium in an external electromagnetic field. He predicted that the contributions of the 
quantum vacuum waves (or "virtual photons") could transfer a nontrivial momentum to 
matter. This prediction was made with the help of the non-relativistic formalism. In our 
opinion, a proper relativistic analysis is needed for a better understanding of the physics 
and of the viability of this phenomenon. Here we begin to reconsider this problem in a 



2 



covariant framework as developed earlier in 4, [ill]. As a first step, we develop a variational 
approach to the description of the dynamics of a moving magnetoelectric medium. The 
corresponding energy-momentum of matter plus electromagnetic field that arises can be 
derived straightforwardly in this formalism from the variation of the total action with respect 
to the spacetime metric. 



II. PRELIMINARIES: THE ESSENCE OF THE FEIGEL EFFECT 

The Feigel effect 



13| can be described in simple terms as follows: Let us consider an 
isotropic homogeneous medium with the electric and magnetic constants e, /i. Electromag- 
netic waves are propagating in such a medium absolutely symmetrically, with the Fresnel 
equation describing the unique light cone. This is easily derived from the constitutive rela- 
tions V = eeq E and Ti = {fi^o)^^B. 

However, if a medium is placed in crossed constant external electric and magnetic 
fields, then it acquires magnetoelectric properties. As a result, we have the anisotropic 
magnetoelectric medium with e, /i, plus the magnetoelectric matrix f] (determined by 
the external fields) which modifies the constitutive relations to V = eeq E + (3 ■ B and 
Ti = {fifio)^^B — ■ E; here denotes the transposed matrix. 

Accordingly, the wave propagation in such a medium also becomes anisotropic and bi- 
refringent, with the wave covectors now belonging to two light cones. Applying this to 
vacuum waves (or, perhaps, better to say to the "vacuum fiuctuations" or "virtual photons" ) 
propagating in the mag netoeleetnc bod. Keigel Q computed the .ota. .no.eMu. ea„ied 
by these waves and concluded that it is nontrivial. In accordance with this derivation, a 



body should move with a sma. 



discussed in 



Q,[l5|, 



16 



17 



18|. 



1 but non-negligeable velocity. Earlier the Feigel process was 



In order to evaluate the possible Feigel effect, it is necessary to substitute the "vacuum 
waves" into the energy-momentum tensor. This paper is devoted to the derivation of the 
latter in the framework of a variational approach. 
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III. CONSTITUTIVE RELATION 



Within the axiomatics of the premetric generally covariant framework the projection 
technique is used to define the electric and magnetic phenomena in an arbitrarily moving 



medium. As in j4|, we assume that the spacetime is foliated into spatial slices with time a 
and transverse vector field n. 

When applying the projection technique to the 2-forms of the electromagnetic excitation 
H and and the electromagnetic field strength F, we obtain the three-dimensional objects: 
the magnetic and electric excitations 7i and V as longitudinal and transversal parts of H 
and, similarly, electric and magnetic fields E and B as longitudinal and transversal parts of 
F, respectively, namely 

H = -nAda + V and F = EAda + B. (3.1) 

This foliation is called the laboratory foliation, with the coordinate time variable a labeling 
the slices of this foliation. 

The spacetime metric g introduces the scalar product in the tangent space and defines 
the line element. With respect to the laboratory foliation coframe it reads (a, b, ... = 1, 2, 3) 

ds"^ = da"^ + gab dx^ = N"^ da^ - gab dx"" dx\ (3.2) 

Here A^^ = g(n, n) is the length square of the foliation vector field n, and daf = dx^—n"' da is 
the transversal 3-covector basis, in accordance with the definitions above. The 3-metric ^^^gab 
is the positive definite Riemannian metric on the spatial 3-dimensional slices corresponding 
to fixed values of the time a. This metric defines the 3-dimensional Hodge duality operator 

The constitutive relation which links the electromagnetic field strength to the electro- 
magnetic excitation, H = H{F), can be nonlocal and nonlinear, in general. Here we will 
confine our attention to the local and linear constitutive relation. 

Then, if we write the the excitation 2-form in terms of its components in a local coordinate 
system {a;*}, (H, V) = H = Hij dx* A dx^ /2 (with = 0, 1, 2, 3), the local and linear 

constitutive relation means that the components of the excitation are local linear functions 
of the components of the field strength {E, B) = F = Fij dx'^ A dx^ /2: 

H = k{F), H,, = ^K,,^'F,i. (3.3) 



Along with the original constitutive K-tensor, it is convenient to introduce an alternative 
representation of the constitutive tensor: 

x'^'':=^e*^""K„„'=^ (3.4) 

Performing a (1 + 3)-decomposition of covariant electrodynamics, as described above, we 
can write H and F as column 6-vectors with the components built from the magnetic and 
electric excitation 3-vectors Tiai'D"' and the electric and magnetic field strengths Ea,B"', 
respectively. Then the linear spacetime relation (13.31) reads: 

(3.5) 

Here the constitutive tensor is conveniently represented by the 6 x 6-matrix 

-i^= r': , = I - - I . (3.6) 







Assuming that the skewon and the axion pieces are absent, we find that the constitutive 
matrices satisfy A"'' = A^", Bab = Bba, and V^" = with Ca = 0. 

The dynamics of a material medium is encoded in the structure of another foliation (r, u) 
which is determined by the four-vector field of the velocity u of matter and the proper time 
coordinate r. Accordingly, we have to formulate the constitutive law with respect to this, 
so called material foliation. As a first step, we observe that the relation between the two 



coframe bases, namely t 
{dr^dx"") is as follows {4]: 



lose of the laboratory foliation {da, daf) and of the material foliation 




-fc/N 


Vb/{cN) 


7t>" 






(3.7) 



Here, for the relative velocity 3-vector, we introduced the notation 



Substituting (13. 7p into (13.21) . we find for the line element in terms of the new variables 



ds^ = dr"^ - gab dx" dx^, where gab = ^^''Qab ^ VaVb- (3.9) 



c 
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The metric ^ab of the material fohation has the inverse 

2 

= (3)^«f' + 2_yV, (3.10) 

where 

(3) gab jg 

the inverse of ^^'Qat- For the determinant one finds (det^ab) = (det ^^'gat) 7 ^• 
We will denote the 3-dimensional Hodge star defined by the metric as *. 

Given the transformation ( 13 .7^ between the 1-form bases of the two foliations (laboratory 
and material), it is straightforward to calculate the components of the constitutive tensor 
with respect to the moving reference frame from the constitutive tensor which describes 
the matter at rest. The corresponding explicit general transformation formulas for the 3- 
dimensional matrices A, B, C can be found in [4]. A nontrivial outcome of this formalism is 
that the isotropic moving medium is described by the constitutive relation based on the so- 
called optical metric (first introduced by Gordon [19|). We will use this fact in the subsequent 
derivations. 



IV. RELATIVISTIC FLUID 



The earlier work includes that for the ideal fluids 



electromagnetic fields was discussed in 23 



24 



m 



27 



Jo 



21 



22| and the case of the fluids in 



26| . Recently, a new study was reported 



An ideal fluid which consists of structure-less elements (particles) is characterized in the 
Eulerian approach by the fluid 4- velocity m*, the internal energy density p, the particle 
density u, the entropy density s, and the identity (Lin) coordinate X. Normally, it is 
assumed that the motion of a fluid is such that the number of particles is constant and that 
the entropy and the identity of the elements is conserved. In other words. 



Vi{uu') = 0, 
u'diS = 0, 



u'diX = 0. 



(4.1) 
(4.2) 
(4.3) 



By the conservation of the entropy only the reversible processes are allowed. In a variational 
approach to continuous media, these assumptions are considered as constraints imposed on 
the dynamics of the fluid by means of Lagrange multipliers. Then the fluid Lagrangian 
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20l . l2l[ | reads Vaiat = -f'mat V--g da A'e with (? = spatial volume 3-form) 



i^mat = - + Ao(t^'^i - c') + Ai V.(Z/M^) + A2 u'diS + A3 m'S^X. (4.4) 

c 

The Lagrange multipliers Ai,A2,A3 impose the constraints (14. II) - (14 .Sp on the dynamics of 
the fluid, whereas Aq provides the standard normalization condition for the 4-velocity 

gijuV = c\ (4.5) 

For the description of the thermodynamical properties of the fluid, the usual thermodynam- 
ical law ( "Gibbs relation" ) is used, 

Tds = d{p/u)+pd{l/u), (4.6) 

where T is the temperature and p the pressure. 

When the medium possesses dielectric and magnetic properties, then we have to add the 
Lagrangian of the electromagnetic field 

Km = - ^H AF = - ^da AiHAB -V AE) 

= -^daAe {A'^'E^E, + Ba^B'^B' - 2C\EaB') . (4.7) 

At first, let us consider the moving isotropic dielectric fluid without magnetoelectric 
properties, that is, C^b = 0. Then the electromagnetic field Lagrangian reads Vem = 
Lem \f--9 da A't with 



1p 

Here Aq is the vacuum admittance, and 



-^^em - - 77 fi'oit 5'opt FikFji. (4.8) 



9l. = 9'' —u'u^ (4.9) 



is the optical metric of Gordon [191]. 



A. Equations of motion without a magnetoelectric effect 

At first, we consider the case when the fluid does not have magnetoelectric properties: 
C"{, = 0. We also assume that there are no free charges and currents in the fluid. 

The dynamics of the system (fluid+field) is thus described by the total Lagrangian V = 
Vmat + Km- The Corresponding equations of motion are then derived from the variational 
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principle for V with the independent variables u\ z/, s, X, Aq, Ai, A2, A3, and Ai (the covector 
of the electromagnetic potential). 

Variation with respect to the Lagrange multipliers Ao,Ai,A2,A3 yields the constraints 
fl4.ip - fl4.3l) and (14.51) . whereas variation of V with respect to u, s,X yields, respectively, 

uu'diAi + = 0, (4.10) 

c 

rri 

V. {A2U') + — = 0, (4.11) 

c 

Vi{Asu') = 0. (4.12) 

In order to derive these results, we used the thermodynamical law (14.60 . from which one has 
dp/ds = vT and dp/du = {p + p)/^. 

Finally, the variation of V with respect to the 4-velocity yields 

2AoUi - udiAi + A^d.s + A^d.X + F^, = 0. (4.13) 

jJiC 

The last term emerges from Vem- Contracting the last equation with and making use of 
dUD-dlJD, dS]), and fICTD . we find 



2A„ = - M^LtP + Mi^ F„uVu, ]. (4. 14) 



The variation with respect to gives the Maxwell equations for the moving fluid 

yA9ilt9%Fki)=Q. (4.15) 



When the fluid has free currents and charges, we have to add the interaction term AjJ* to 
the electromagnetic field Lagrangian Lem, and then the right-hand side of (14.151) picks up a 
nontrivial current J*. 



B. Energy-momentum tensor 

The above consideration is not only relativistically covariant, but actually generally co- 
variant. The metric gij above describes an arbitrary curved spacetime. The invariance of the 
total Lagrangian L under general coordinate transformations (using the standard Noether 
machinery) yields the conservation of the energy-momentum of the system (fiuid+field) . We 
obtain the energy-momentum tensor as usual from the variation of the total Lagrangian with 
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respect to the metric, i.e., 



Tij . 2 c 



mat ~r J-'e 



The computation is straightforward, and it yields for (14.41) and (14. 8p 



(4.16) 



1 



p + p 



UiU-' + 



+ 



k TjrmL 



(4.17) 



In the course of its derivation, we made use of the resuhs of the previous subsection, in 
particular of the equations (I4.10p and (14.140 . The first line in (14.170 originates from Knat, 
whereas the rest comes from Vem- It is interesting though that it is possible to rearrange the 
contributions of the fluid and of the field in such a way that the total energy-momentum is 
recast into the nice form 



Pes + 



Pes + PeS 



■ UiU-' 



+ 



1 



1 • 1 — r?^ 



4 - c- 

Here we introduced the effective pressure and the effective energy density by means of 

l-n2 



(4.18) 



Pes 



P + 



Pes ■= P 



2ppo c2 



Fiu'F'W, 
F^u^F^'uk. 



(4.19) 
(4.20) 



2ppo 

General covariance underlies the conservation law of the total energy-momentum tensor 



V,T,^ =0, 



(4.21) 



which yields the equations of motion of the medium. 



C. Extension to the magnetoelectric case 

The above results all refer to the isotropic dielectric and magnetic medium in motion 
without magnetoelectric crossterms. The extension to the magnetoelectric case, see O'Dell 
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291], is straightforward. The electromagnetic field Lagrangian (14.80 is replaced by the more 



general one 

^em = - ^ 4t ^opt F^kFji - ^ (3'' J, F,,U' V^'^'^^FimUk. (4.22) 

Here the tensor (3^j gives the relativistic generalization of the magnetoelectric matrix (3. The 
totally antisymmetric Levi-Civita tensor is defined as usual by 77'-''^' = e*-''^'/ \^—g, with the 
permutation symbol e^^^'' that have the only nontrivial component, e°^^^ = 1. The dimension 
[P'j] = [Ao]. 

The tensor is assumed to be traceless. If a nontrivial trace is included, this will describe 
an axion contribution. Indeed, the trace = 6j a, when substituted into fl4.22p . yields the 
term ~ aFinu^rf'^^^^FimUk- Using the identity (which states that a totally antisymmetric 
tensor of the 5-th rank is zero in four dimensions) 2u^'^7]^'^^^'^ = u'^r]^'^'-'^"- + y}r]'^'^'^^ + ^i™?^*'^'*^ 
we find a FiiiU^rf'^^'^FimUk = ^arf^^^FijFu- Thus, indeed the trace 5 adds an axion con- 
tribution to the electromagnetic Lagrangian. We assume that the axion is absent in this 
paper. 

Furthermore, it is easy to see that only a projection of the magnetoelectric matrix on the 
rest frame of the 4- velocity u enters the Lagrangian (14.221) . namely 

- ^ F,,v}r^^'^''^FimU, = ^''y F.iu'r^^'^'^Fi^Uk^ (4.23) 



where 



(3'^ = PIP P'u (4.24) 



with the projector defined by PI := 61 — u^Uk/c^- By construction, u^P^j = and UiP'^j = 0, 
which means that only the 3-dimensional transversal part of the magnetoelectric tensor 
contributes. 

Computation of the electromagnetic excitation H = — (9Km/ dF is straightforward. In the 
(l-|-3)-decomposed form, this yields the constitutive relation (13. 5p . where the 3-dimensional 
matrices read explicitly: 

^ab ^ _lab_ ^^ab^ (4 25) 

o -, me -, 

Bab = (/^-')afe + (/^-')afe, (4.26) 
o me 

C\ = Tb+ 7V (4.27) 



Here the first terms on the right-hand sides describe the contributions of the isotropic moving 
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medium 



9 



O 1 



c(l — f ^/c^) 
Ac 



n 



ab 



(3) 



1 nv^ 



1 

72 

n 



(4.28) 
(4.29) 
(4.30) 



1 — f 

Here A := ^Jeeq/ and n := y^e/J is the refraction index of the medium. The 3-dimensional 
indices are raised and lowered with the help of the 3-dimensional spacetime metric ^^^gab, 
see the line element (13. 9p . 

The magnetoelectric contributions to the constitutive matrices read 

2 



/i- 



1 — v'^/c^ — 
2 



7' 



Here we introduced the notation 



1 — v'^/c^ 



(^{a ^b)cdV'^i 



1 — f ^/c^ 



(4.31) 
(4.32) 
(4.33) 



(4.34) 



(4.35) 



with the "doubly projected" magnetoelectric matrix 

For the medium at rest (with f = 0), we straightforwardly verify that the above formulas 
reduce to the constitutive law 



Ha = — (3^a Ef, H Ba- 



(4.36) 
(4.37) 



From these, the constitutive relations of the moving magnetoelectric medium fl4.25l) - fl4.33p 
can be alternatively derived with the help of the transformation (13. 7p that links the labora- 
tory to the material foliations. 

The new magnetoelectric term in the Lagrangian fl4.22p will modify the equations of 
motion of the medium and the total energy-momentum of the system. Direct computation 
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yields: 



rr j d I ^eff + Pes j , 2 x 

Ti^ = - Oi Pes H 5 UiU^ + - quu-" 

& c 



1 

+ — 



2 



1 1 ri 

4 



(4.38) 

yUo/U [ 4 " J 

Here the effective energy density and the effective pressure remain the same fl4.19p . fl4.20p . 
whereas the vector is defined by 

(li = \ P'f PH'I'^Fi^F.kuK (4.39) 

Clearly, it is orthogonal to the velocity of the medium, u^qi = 0. The terms with the 
structure 2u(iqj)/c are well known in the models of relativistic fluids; they usually describe 
fluids with fluxes of energy. It is interesting that the magnetoelectric parameters introduce 
such an additional energy flux term in the total energy-momentum tensor. 

At the beginning, we assumed that the axion is absent. However, it is instructive to 
check how it could contribute to the energy-momentum. For this, we relax for a moment the 
tracefree condition for the magnetoelectric tensor and allow for a nontrivial trace = 5* 5. 
This yields a contribution ~ ^aP-rj^'/'^FimFjikU^ in (14.391) . With the help of the identity 
2u^^r]^'^j-^ = Uj'T]^'^^'^ + u^'r]^' + u'^r]^' , this is transformed into P^Ujr]^'^''"^FimFjik- It 
vanishes because of the identity -P/wj = 0. In other worlds, this demonstrates that the axion 
does not contribute to the energy-momentum, in accordance with its general theory jj]. 



D. Explicit components of the energy-momentum 

All the results of the previous two subsections are generally covariant, i.e., they are 
valid in an arbitrary curved spacetime with any metric gij. Now, returning to a more 
narrow case actually discussed by Feigel, we will specialize to the Minkowski spacetime with 
gij = diag(c^ -1, -1, -1). 

In order to write down the separate components of the energy-momentum tensor, we first 
notice that = 'y,u"' = 7f", with the 3- velocity and the usual relativistic (Lorentz) 
factor 7 = 1/ a/1 — v'^/c^. Then, denoting Xi := FijU^ , we find 

Xo = -i{vE), X, = 7(E + vxB),. (4.40) 

Accordingly, we have the invariant [that enters (I4.19p . (I4.20p 

XiX' = -7^ [E^ - {vEf/c^ + v'^B'^ - {vBf + 2(E ■ [v x B])] . (4.41) 
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Here we denoted the scalar product (vE) = (v ■ E) and similarly (vB). 
After these preliminaries, we find, componentwise, the energy density 



rp 

J-0 



M = 7 



Pcflf + ^Peflf + 



2/1 \ Ho 



[vEf 



9a, 



(4.42) 



the energy flux density (or Poynting vector) 



1q — S 



7 



(Peff + Peff) V + (^i^)(E + V X B) 



+ ^ [E X B]'^ - 7C (st - -yv^ q\ (4.43) 



the momentum density 

-'^ a Fa 



1 



V + ^^^il^(.E)(E + V X B) 



^ [B X E]„ + ^ - \v'v^ q,, (4.44) 
fl c \ c / 



and the stress tensor 



I Xb , 'y b \ 2Pett b 

-Pea + ^VaV - 7 — VaV 



+ 



/i V 2 " 

£0(1 - 



iPeff 

c' 

1 



1 



fiflQ \ 2 



■ 7^ (E + V X B),(E + V X B)^ - 2^q^a' 



(4.45) 



As we see, the vector indeed induces an additional energy flux in (14.431) and adds a 
corresponding magnetoelectric term in the electromagnetic momentum density (14.441) . 

For completeness, let us also give the explicit expressions for the effective energy density 
and the effective pressure of the medium: 



Peflf = P 
PcS = P 



£0(1 -ri2) 2 



2^ 



-7 



2/i 



7 



E^ - {vEf/c^ + v^B^ - {vBf + 2(E • [v x B]) , (4.46) 
E^ - {vEf/c^ + v'^B'' - {vBf + 2(E ■ [v x B])l . (4.47) 



The energy-momentum tensor obtained here is different from both the Minkowski and 
the Abraham energy-momentum tensors. It is thus necessary to check whether the new 
expression is consistent with the main experiments in phenomenological electrodynamics. 
As a first step, we consider usual matter without magnetoelectric properties. 
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V. EXPERIMENTS WITH LIGHT PRESSURE 



In order to test how the formal 
iments of Richards and Jones 30 , 



ism works, we analyse in this section the famous exper- 



3l| . Their measurements have shown that the radiation 
pressure on a metalhc plate immersed into a dielectric fluid is proportional to the refraction 
index n of the medium. 

The problem of the propagation of a plane wave in a dielectric medium is solved exactly. 
Suppose a wave travels along the x axis, with the metal plate surface located at x = 0. 
For a; < 0, we have dielectric fluid with the permeability and permittivity and for 

X > we have the metal characterized by 82-, /i2 and by the electric conductivity a. Then the 
electric and magnetic fields that describe the plane waves normally incident on and reflected 
form the metal surface are (x < 0) 



a e 



-iijjt 



iDz = — CL e 



+ 7^e 



The fields in the metal (x > 0) read 



n2 
c 



1 



(5.1) 
(5.2) 

(5.3) 
(5.4) 



Without loosing generality, we assume that the electric field is directed along the y axis. 
We use the complex representation to simplify the formulas. Here the complex variable a 
describes the amplitude of the incident wave, whereas the complex reflection and transmis- 
sion coefficients are denoted 71 and T, respectively. As usual, rii = ^/sipi and n2 = ^£2/^2- 
Solving the Maxwell equations, we find IC = k + i6 with 






1/2 










)] 











k 



From the matching conditions at the boundary x = 0, we obtain explicitly 




1/2 



(5.5) 



n 



/c 



r 



2, 

Ml 



/c 



(5.6) 



All the formulas are simplified in the case of the strongly conducting matter, when a ^ ssq uj 
(see [lo|, e.g.), but this condition is actually not necessary in our analysis. 
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It is now straightforward to compute the force per area of the metal surface, which is 
most conveniently given in terms of the stress tensor components, fx = — Sx^ ■ We put 
f ° = since the dielectric fluid is at rest. Then fl4.45p yields 

2 2/ii/io 
= P+^(l + |7^P).^,). (5.7) 

Here the right-hand side is evaluated at x = 0, and we used (14.431) to substitute the energy 
flux (Poynting vector) of the incident wave s^-^ along the x axis. As usual, the averaging 
over a period was performed in these calculations. The hydrodynamical pressure term p 
is obviously cancelled by the same fluid pressure force acting on the metal plate from the 
opposite side. So, the net force per area on the is described only by the last term in (I5.7j) . in 



complete agreement with the experimental results (cf. also with the formula (4.22) in 10|). 



Another application o 
Ashkin and Dziedzic 32, 



our formalism to the analysis of the important experiments of 



33| will be considered elsewhere. 



VI. DISCUSSION AND CONCLUSION 

In this paper we have developed a model for a moving electrodynamical medium within 
the framework of a variational approach. The application to the observations of Richards 



and Jones 



30 



31 1 demonstrates the viability of the model. Our next aim is to use this 

n 

model for the analysis of the possibility of the Feigel effect pj| predicted recently. 

We have derived the total energy-momentum of moving matter and the electromagnetic 
field. In particular, we found that the magnetoelectric properties contribute with the vector 
Qi to an additional energy flux and momentum density. 

We are now in a position to compare the quantities derived here with those used by 



Feigel [13[ . A direct inspection shows that the expression of Feigel is not obtained from the 
formulas above in the non-relativistic limit, when we assume that v"^ /(? <^ 1. This fact casts 
certain doubts on the validity of the original conclusions of Feigel. 



In the meantime, an interesting new paper [3J| appeared where the Feigel effect was 
analysed with the help of a Green-function techniques. A similar discussion, which is based 
on the quantum-mechanical formalism and using the above variational formalism, will be 
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published elsewhere. 
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